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ABSTRACT 



fV^ Context, We study the circumstellar environment of the carbon-rich star R Sculptoris using the near- and mid-infrared high spatial 

ry^ \ resolution observations from the ESO-VLTI focal instruments VINCI and MIDI, respectively. 

Aims. These observations aim at increasing our knowledge of the dynamic processes in play within the very close circumstellar envi- 
i-Q ■ ronment where the mass loss of AGB stars is initiated. 

Methods. We first compare the spectro-interferometric measurements of the star at different epochs to detect the dynamic signatures 
of the circumstellar structures at different spatial and spectral scales. We then interpret these data using a self-consistent dynamic 
model atmosphere to discuss the dynamic picture deduced from the observations. Since the hydrodynamic computation needs stellar 
parameters as input, a considerable effort is first applied to determining these parameters. 

Results. Interferometric observations do not show any significant variability effect at the 16 m baseline between phases 0.17 and 0.23 
in the K band, and for both the 15 m baseline between phases 0.66 and 0.97 and the 31m baseline between phases 0.90 and 0.97 in the 
N band. We find fairly good agreement between the dynamic model and the spectrophotometric data from 0.4 to 25 /im. The model 
agrees well with the time-dependent flux data at 8.5 fim, whereas it is too faint at 11.3 and 12.5 /mi. The VINCI visibility measure- 
ments are reproduced well, meaning that the extension of the model is suitable in the K-band. In the mid-infrared, the model has the 
proper extension to reveal molecular structures of C 2 H 2 and HCN located above the stellar photosphere. However, the windless model 
, used is not able to reproduce the more extended and dense dusty environment. 

C^) • Conclusions. Among the different explanations for the discrepancy between the model and the measurements, the strong nonequilib- 

rium process of dust formation is one of the most probable. The transition from windless atmospheres to models with considerable 
mass-loss rates occurs in a very narrow range of stellar parameters, especially for the effective temperature, the C/O ratio, and the 
pulsation amplitude. A denser sampling of such critical regions of the parameter space with additional models might lead to a better 
representation of the extended structures of low mass-loss carbon stars like R Sculptoris. The complete dynamic coupling of gas and 
dust and the approximation of grain opacities with the small-particle limit in the dynamic calculation could also contribute to the 
difference between the model and the data. 



Key words, techniques: interferometric - techniques: high angular resolution - stars: AGB and post-AGB - stars: atmosphere - stars: 
circumstellar matter - stars: mass-loss 

1 . Introduction As the opacity of amorphous carbon dust is high, grains receive 

enough momentum through radiative pressure to be accelerated 

AGB stars are complex objects where several physical processes and tQ dmg along ±e gafj by collisionSi causing a slow outflow 

are in play (nucleosynthesis, convective effects, stellar pulsation, from the star ( F i e i sc her et al. 1992; Hofner & Dorfi 1997). 
shock waves, mass-loss, molecular and dust formation, etc.). To investigate this dynamic process , it is necessary to com- 
Among these mechanisms, the mass loss is particularly impor- pare predictions of hydrodynamic atmospheric models with 
tant for the AGB phase by limiting the maximum luminosity observa tions. The code described by Hofner et al. (2003) 
of the object and controlling the overall stellar evolution. The solyes the coupled equations of hydrodynamics, together 
general picture of the mass-loss process of carbon-rich AGB wkh frequency-dependent radiative transfer, taking the time- 
stars is quite well understood nowadays. The pulsation of the dependent formation; growth , an d evaporation of dust grains 
star creates shock waves pushing the matter to conditions where intQ account Dynamic mode i atmospheres coming from this 
both temperature and pressure allow formation of dust grains. code have successfully re p rod uced line profile variations 

(Nowotny et al., 2005a,b) and time-dependent spectrometric 

Send offprint requests to: S. Sacuto data (Gautsch y-Loidl et al., 2004) of carbon-rich stars. 

* Based on observations made with the Very Large Telescope Because of the large distance of the closest AGB stars (typically 

Interferometer at Paranal Observatory under programs 60.A-9220, lnn , mnn \ i.uu-u i w .. u ■ 

074.D-0601, 077.D-0294 (French Guaranteed T.me Observation), 1 9°, to 10 u ° P al " SeC) ' ^ * e hl f ^ Solutl0n ^ chm ^ e 

078.D-0112 (Belgian Guaranteed Time Observation), and 078.D-0122 of long-baseline optical interferometry allows probing the 

(French Guaranteed Time Observation). regions where the mass-loss process develops. Interferometric 
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predictions based on the dynamic atmospheric models with 
and without mass loss are described in detail by Paladini et al. 
(2009). The aim of this paper is to extend the work comparing 
the dynamic models with near- and mid-infrared high angular 
resolution interferometric observations of a carbon-rich star. 
RSculptoris (RScl) is a carbon-rich, semi-regular (SRa) vari- 
able star having a low mass-loss rate (l-5xl0~ 7 M yr _1 : Le 
Bertre 1997; Gustafsson et al. 1997; Wong et al. 2004). The 
mid-infrared emission between 10 and 12 //m observed in the 
ISO/SWS data indicates that the star is surrounded by a warm, 
carbon-rich dusty shell composed of amorphous carbon (AmC) 
and silicon carbide (SiC) (Hron et al., 1998). The pulsation 
period of the star is 374 days (Whitelock et al., 1997), and the 
cycle-to-cycle averaged visual magnitudes vary from 6.7 to 8.1. 
A second period of 1804 days (~5 times the stellar pulsation 
period) has also been detected by Whitelock et al. (1997); 
however, its origin is not yet understood. 

Until now, all the interferometric studies of temporal evolution 
of carbon-rich stars have been based on sequences of geo- 
metric and hydrostatic models with varying stellar parameters 
(van Belle et al., 1997; Thompson et al., 2002; Ohnaka et al., 
2007). However, as argued by Hofner et al. (2003), the only 
way to understand how the dynamic processes influence the 
atmospheric structure at different spatial scales is the use of 
time-dependent, self-consistent dynamic models. 
In this paper, we thus present the first interpretation of spec- 
trophotometric and interferometric data of a carbon-rich star in 
terms of a self-consistent dynamic model atmosphere. Near- and 
mid-infrared interferometric data at different baselines (15, 30, 
and 60 m) will allow probes of regions ranging from the stellar 
photosphere to the innermost region of its dusty environment 
where the dynamic processes of shock waves, stellar winds, and 
mass loss are all at work. 

The outline of this paper is as follows. In Sect. 2, we present 
the VINCI and MIDI observations of R Scl and describe the cal- 
ibration performed with the corresponding data-reduction soft- 
ware packages developed for these instruments. Section 3 is 
dedicated to the variability of the star obtained from previous 
visual photometric and polarimetric data. We investigate the 
mid-infrared spectrometric and near- and mid-infrared interfer- 
ometric variability of the star by comparing measurements from 
cycle-to-cycle and phase-to-phase. In Sect. 4, we investigate the 
hydrodynamic modeling of R Scl. Since the hydrodynamic com- 
putation needs stellar parameters as input, a considerable effort 
is first applied to the determination of the corresponding stellar 
parameters. The last part of this section presents the results of 
the self-consistent hydrodynamic modeling of the spectrophoto- 
metric and of near- and mid-infrared interferometric data of the 
star. Finally, we conclude and give perspectives for future work 
in Sect. 5. 

2. Observations 

The Very Large Telescope Interferometer (VLTI) of ESO's 
Paranal Observatory was used with VINCI, the near-infrared 
(/l=2.0-2.4 fim) interferometric recombiner (Kervella et al., 
2000), and MIDI, the mid-infrared (/l=8.0-13 /mi) interferomet- 
ric recombiner (Leinert et al., 2003). VINCI operates around 
2 microns (K band), using single spatial mode interference by 
employing optical fibers for spatial filtering and beam combina- 
tion. MIDI combines the light of two telescopes and provides 
spectrally resolved visibilities in the N band atmospheric 
window. 



Observations of R Scl were conducted in 2001, 2005, and 2006 
with the VLT unit telescopes (UTs) UT2-UT4, and auxiliary 
telescopes (ATs) E0-G0, D0-G0, A0-G0, G0-K0, G0-H0, and 
D0-H0. This provides projected baselines and projected angles 
in the range of 11 to 64 m and 34° to 117°, respectively (see 
Tables 1, 2, and 3, and Fig. 1). The seeing of all the VINCI 
observations is below 0.9", whereas it is below 2.4" in the case 
of the MIDI observations (see Tables 2 and 3). 
Tables 1, 2, and 3 present the journals of interferometric 
observations of RScl with the VINCI and MIDI focal in- 
struments. The calibrators, jSCet (G9II-III, angular diameter: 
0=5.18+0.06 mas), aCMa (A1V, 0=5.94±O.O2 mas) for 
VINCI and HD12524 (K5III, 0=2.53±O.O7 mas), /SGru 
(M5III, 0=26.8+1.3 mas), T 4 Eri (M3/M4III, 0=10.6+1.0 mas), 
HD4128 (G9II-IH, 0=5.04+0.07 mas) for MIDI, were observed 
before or after the science target. Diameters of the VINCI 
calibrators were adopted from the CHARM2 catalog (Richichi 
et al., 2005). Diameters of the MIDI calibrators HD12524 and 
HD4128 were taken from the ESO website through the CalVin 
database 1 . The diameter of the MIDI calibrator t 4 Eri was taken 
from the CHARM2 catalog (Richichi et al., 2005). Finally, the 
diameter of the MIDI calibrator /3Gm was defined from the 
work of Sacuto et al. (2008). 

The VINCI data (program 60.A-9220) were processed with 
the pipeline of data reduction based on the wavelet transform 
described in Kervella et al. (2004). After checking the stability 
of the transfer function (r,) over the two nights of observations, 
the raw visibility measurements were calibrated using the 
weighted average of the estimated Tf. 

The VINCI instrument has a bandpass corresponding to the 
K-band filter (2.0-2.4 yum) without spectral dispersion. It is thus 
necessary to compute an effective wavelength to determine the 
spatial frequency of the observation (Kervella et al., 2007). 
The effective wavelength of our observations was calculated 
from the relation given by Wittkowski et al. (2004) using the 
best-fitting hydrostatic model spectrum of RScl described in 
Sect. 4.1. We obtained an effective wavelength of 2.253 //m. 
Among the VINCI data, one observation (2001-12-08 01:03:54; 
see Table 1) led to an outlier showing an abnormally low level 
of amplitude as compared to the other visibility data at close 
base, PA, and phase. The corresponding visibility measurement 
was discarded. 

In the case of MIDI, chopped acquisition images were recorded 
(f=2Hz, 2000 frames, 4 ms per frame) to ensure the accurate 
acquisition of the target. Photometry from the UT observations 
(program 074.D-0601) was obtained before and after the inter- 
ferometric observations with the HIGH-SENS mode of MIDI, 
using the GRISM that provides a spectral resolution of about 
230. Photometry from the AT observations were obtained with 
the HIGH-SENS mode (program 078.D-0112), and simultane- 
ously with the interferometric observations with the SCI-PHOT 
mode of MIDI (programs 077.D-294 and 078.D-0122), using 
the PRISM that provides a spectral resolution of about 30. 
The acquisition images enabled us to confirm that the mid-IR 
source was unresolved by the single-dish UT. This implies that 
most of the mid-infrared flux originates in the inner 300 mas 
of the source, corresponding to the Airy pattern of the UTs 
and defining the Field Of View (FOV) of the interferometric 
observations. This also means that the flux detected within the 
1.1" FOV of the single-dish AT originates in this inner region. 
Therefore, measurements extracted from the ATs and the UTs 



1 http : //www . eso . org/observing/etc/ 
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# 


Star 


UT date & Time 


Phase 


Config. 


Basefm] 


PA[deg] 


Seeingfarcsec] 


Mode 


Spectral resolution 


1 


RScl 


2005-01-03 03:01:23 


0.23 


U2-U4 


61 


117 


1.7 


HIGH-SENS 


230 


Cal 


HD12524 


2005-01-03 03:27:24 










1.0 






2 


RScl T 


2005-01-03 04:05:54 


0.23 


U2-U4 


47 


136 


0.9 


HIGH-SENS 


230 


Cal 


HD 12524 


2005-01-03 04:45:32 










0.8 







Table 3. Journal of all available MIDI Auxiliary Telescopes observations of R Scl. 



# 


Star 


TTT Hatp # Timp 


Phase 


Confi 


Rocprml 
ud&CLinj 


PA TrWl 


\ppin cj TaTr cpr 1 
occiiigLciiuftCLj 


Mode 


Sv\prtTtn rpcr\ln li r\n 

OlJCLlldl 1CMJ1UL1U11 


3 


RScl 


2006-06-16 09 38 09 


0.64 


E0-G0 


14 


42 


0.7 


SCI-PHOT 


30 




R dm 
P KJl 11 


900ft oft \f\ 10-91-fto 










O ft 
u.u 






4 


R Scl 


2006-06-7? 09-29-58 


0.66 


D0-G0 


29 


46 


0.8 


SCI-PHOT 


30 


Cal 




900ft-0ft-99 09-49-48 










0.9 






5 


R Scl 1 ' 


2006-06-23 09:23:50 


0.66 


E0-G0 


14 


45 


1.0 


SCI-PHOT 


30 


5 


R Scl 


7006-06-71 09-46 04 


0.66 


E0-G0 


15 


50 


1.4 


v > v i i n w i 


30 


Cal 


T 4 Eri 


2006-06-23 10:04:33 










1.1 






7 


R Scl 


7006-08-08 0Q'57'10 


0.79 


A0-G0 


62 


77 


1.1 


iH i i nw i 


30 


Cal 


T 4 Eri 


2006-08-08 10:09:22 










1.4 






Cal 


nut izo 


900ft OQ 17 04-5ft-51 










9 O 
z.u 






8 


R Scl "'" 


2006-09-17 05:26:12 


0.89 


A0-G0 


63 


63 


1.8 


HIGH-SENS 


30 


i^ai 


iLUH 1 Zo 


900ft OQ 17 0S-S1-4^ 
zuuo-uy-i/ uj.ji.'tj 










1 7 






9 


R Scl 1 " 


2006-09-17 06:20:53 


0.89 


A0-G0 


64 


71 


2.4 


HIGH-SENS 


30 


Til 


Hn/M 9« 

n i zo 


900ft 09 17 Oft-49-41 
ZUUO-U:? 1 / UO.4-Z.4-l 










2 






10 


RScl f 


2006-09-17 07:06:01 


0.89 


A0-G0 


63 


77 


1.9 


HIGH-SENS 


30 




Hn41 98 
n \J-r 1 Zo 


900ft 09 18 O'M'Mft 










1 5 
l . j 






11 


R Scl 1 " 


2006-09-18 03-55-48 


0.89 


A0-G0 


59 


48 


1.6 


HTGH-SFNS 

1 J.XVJ1 1 Lll_il 1 i-t 


30 


12 


R Sri? 
Ix OC1 


900ft 09 1 8 04- 19-1 5 


u.oy 


AO no 


fto 

UU 


si 


1 5 
l . J 


HTHH SFMS 




Cal 


HD4128 


2006-09-18 04:39:39 










1.7 






1 ^ 

1 J 


R Srlt" 
i\ OL1 


900ft OQ 1 Si OS -09 -OS 

zuuo-uy-io uj.uz.uo 


ft SQ 


An no 

nu-uu 


ft9 

oz 


ou 


i ft 

1.0 


HTHH SPMS 


jU 


Cal 


HD4128 


2006-09-18 05:25:49 










1.6 






Til 


Hn41 98 
nut i zo 


900ft 09 18 07-05-14 
ZUUO-U7 - 1 o U/.UJ.1H- 










1 5 






14 


RScl* 


2006-09-18 07-34-30 

—\jyj\j yj j i o yj i . j^t. j\_/ 


0.89 


A0-G0 


61 


80 


2.1 


HTGH-SFNS 

1 J.XVJ1 ± LlJ_il 1 *_J 


30 




Hn41 98 
nu 4 + 1 zo 


900ft 09 19 07 07 0ft 
ZVAfO-U?- 1 7 U/.U/.UO 










1 2 






15 


RScl* 


2006-09-19 07-32-00 


0.90 


D0-G0 


30 


81 


1.3 


HIGH-SENS 


30 


i_ai 


LfD/ll 90 

rlU'+lzo 


900ft OQ 1Q 07-Sft-OI 

zuuo-uy-iy u/.jo.ui 










1 9 
1 .z 






16 


RScl f 


2006-09-19 08:19:20 


0.90 


D0-G0 


28 


87 


1.5 


HIGH-SENS 


30 


Cal 


HD41 98 
n \J-r 1 Zo 


900ft 09 90 04- 19-97 
ZUUO-U7-ZU UM-.IZ.Z/ 










1 3 






17 


RScl 


2006-09-20 04:36:30 


0.90 


D0-G0 


31 


57 


1.6 


HIGH-SENS 


30 


1 R 
I o 


R Sri 
l\ OL1 


900ft 09 91 09-40-17 
z.\)Uyj-Uy- Z 1 UZ.HU. 1 / 


n on 


no KTl 


54 


jt 


1 7 
l . / 


UTf.TJ SPMS 


^0 


Cal 


HD4128 


2006-09-21 03:03:59 










1.1 






1 Q 


R Sri 
I\ OL1 


900ft 10 1ft 0105-9ft 
ZUUU-1U-1U UJ.UJ.ZU 


Q7 
u.y i 


uu-nu 


T1 

j l 


J ^ 


ft 
u.u 


SPT PHOT 
j v i - 1 n v / 1 


^0 


Cal 


Eri 


2006-10-16 03:28:25 










0.5 






ZU 


R Sri 
IN OC1 


900ft 10 1ft 07-00-0^ 
ZUUD-1U-1D U/.UU.Uj 


n Q7 


UU-rlU 


5^ 

j j 


on 


O ft 
U.D 


sr*T PHOT 

jv l-] rlV / 1 


^O 
jU 


Cal 


Eri 


2006-10-16 07:24:05 










0.5 






21 


R Scl 


900ft 10 17 09 05 09 
ZUUO-1U 1 / UZ.UJ.UZ 


97 


D0-H0 


59 


t ^ 


1 4 


SPT PHOT 
. j v i r nv/ 1 


^0 


Cal 




2006-10-17 02:26:29 










1.0 






99 
ZZ 


R Sri 
IN OLl 


900ft 10 18 01-51-95 
ZUUD-lU-lo Ui.jl.Zj 


n Q7 


DU-UU 


1 S 

1 j 


t / 


U.o 


HTHH SPMS 
ri IVi ri -. 1 1 L.N o 


^O 


Cal 


HD4128 


2006-10-18 02:13:41 










0.9 






23 


RScl 


2006-10-18 02:34:58 


0.97 


E0-G0 


15 


55 


0.9 


SCI-PHOT 


30 


Cal 


r 4 Eri 


2006-10-18 02:59:24 










0.8 






24 


RScl 


2006-10-18 04:51:15 


0.97 


E0-G0 


16 


75 


1.0 


HIGH-SENS 


30 


Cal 


HD4128 


2006-10-18 05:22:11 










1.1 






25 


RScl T 


2006-10-19 03:28:15 


0.98 


E0-G0 


16 


64 


1.0 


HIGH-SENS 


30 


Cal 


HD4128 


2006-10-19 03:52:03 










0.4 






26 


RScl 


2006-10-19 04:15:04 


0.98 


E0-G0 


16 


71 


0.7 


HIGH-SENS 


30 


Cal 


HD4128 


2006-10-19 04:37:45 










0.7 






27 


RScl 


2006-10-19 05:23:42 


0.98 


E0-G0 


15 


79 


0.5 


SCI-PHOT 


30 


Cal 


t 4 Eri 


2006-10-19 05:49:48 










0.4 






Cal 


r 4 Eri 


2006-12-16 02:57:16 










1.7 






28 


RScl 


2006-12-16 03:19:24 


0.13 


G0-H0 


25 


93 


1.4 


SCI-PHOT 


30 


Cal 


t 4 Eri 


2006-12-18 00:57:42 










1.2 






29 


RScl 


2006-12-18 01:21:29 


0.14 


G0-H0 


31 


79 


1.2 


SCI-PHOT 


30 


Cal 


r 4 Eri 


2006-12-21 01:27:43 










0.9 






30 


RScl 


2006-12-21 02:06:31 


0.15 


G0-H0 


29 


86 


1.7 


SCI-PHOT 


30 
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Table 1. Journal of all available VINCI observations of R Scl. 



Star 


UT date & Time 


Phase 


Config. 


Base[m] 


PA[deg] 


/jCet 


2001-11-15 01:30:15 




- 






/jCet 


2001-11-15 01:54:23 




- 






ySCet 


2001-11-15 02:01:16 




- 






RScl 


2001-11-15 02:15:27 


0.17 


EO-G0 


16 


69 


RScl 


2001-11-15 02:29:56 


0.17 


E0-G0 


16 


71 


RScl 


2001-11-15 02:35:41 


0.17 


E0-G0 


16 


71 


RScl 


2001-11-15 02:40:42 


0.17 


E0-G0 


16 


72 


/jCet 


2001-11-15 03:50:07 




- 






/jCet 


2001-11-15 03:58:02 




- 






ySCet 


2001-11-15 04:05:08 




- 






RScl 


2001-11-15 04:18:11 


0.17 


E0-G0 


15 


84 


RScl 


2001-11-15 04:34:41 


0.17 


E0-G0 


14 


86 


RScl 


2001-11-15 04:39:59 


0.17 


EO-GO 


14 


87 


/?Cet 


2001-11-15 05:32:55 




- 






/jCet 


2001-11-15 05:39:46 




- 






ySCet 


2001-11-15 05:46:41 




- 






RScl 


2001-11-15 06:03:45 


0.17 


EO-GO 


11 


99 


RScl 


2001-11-15 06:09:24 


0.17 


EO-GO 


11 


100 


RScl 


2001-11-15 06:14:57 


0.17 


EO-GO 


11 


101 


((Cma 


2001-11-15 06:30:06 




- 






ffCma 


2001-11-15 06:34:52 




- 






orCma 


2001-11-15 06:46:09 




- 






ffCma 


2001-11-15 06:52:23 




- 






RScl 


2001-12-08 00:57:37 


0.23 


EO-GO 


16 


70 


RScl f 


2001-12-08 01:03:54 


0.23 


EO-GO 


16 


72 


RScl 


2001-12-08 01:13:17 


0.23 


EO-GO 


16 


72 


RScl 


2001-12-08 01:18:42 


0.23 


EO-GO 


16 


73 


RScl 


2001-12-08 01:25:45 


0.23 


EO-GO 


16 


74 


/jCet 


2001-12-08 01:38:37 










ySCet 


2001-12-08 01:45:31 










/jCet 


2001-12-08 01:50:26 










RScl 


2001-12-08 02:09:03 


0.23 


EO-GO 


15 


79 


RScl 


2001-12-08 02:13:53 


0.23 


EO-GO 


15 


80 


/jCet 


2001-12-08 02:26:18 










/jCet 


2001-12-08 02:33:18 











Note 1 : The calibrators used to calibrate the visibilities are given before 
or after the science target. The visual phase of the star during the 
observations is calculated from Eq. 1. The configuration used for the 
observations is given. The length and position angle of the projected 
baseline are also indicated. Same for Tables 2 and 3. 



Note 2: 1 means that the data were not exploitable (see text). Same for 
Tables 2 and 3. 



are fully consistent with each other. 

The data reduction software packages 2 MIA and EWS (Jaffe, 
2004) were used to derive the calibrated spectra and visibilities 
(Chesneau et al., 2005; Ratzka et al., 2007). MIA is based on a 
power spectrum analysis and uses a fast Fourier transformation 
(FFT) to calculate the Fourier amplitude of the fringe packets, 
while EWS uses a shift-and-add algorithm in the complex 
plane, averaging appropriately modified individual exposures 
(dispersed channeled spectra) to obtain the complex visibility. 
To calibrate the MIDI visibility measurements of RScl, we 
used the transfer functions derived from the closest calibrators 
in time from the source (see Tables 2 and 3). This allows a 
decrease in the bias due to the air mass difference between the 
source observation at a given sky location and the calibrator 
observation at another location. The main source of uncertain- 
ties on the calibrated visibility measurements comes from the 
time-dependent atmospheric transmission between the obser- 
vation of the calibrator and the source and from the dimension 



and uncertainty on the calibrator angular diameter (Borde et al., 
2002; Cruzalebes et al., 2010). The rms of the transfer function 
deduced from all the calibrators of the night gives an upper esti- 
mate of the calibrated visibility error bar due the time-dependent 
atmospheric condition between the observations of the source 
and the calibrator. We assume a standard error bar of 10% on 
the calibrated visibilities corresponding to the upper limit on 
the error budget previously assigned. However, in the case of 
the HIGH-SENS mode, the photometry is performed about 5 
to 10 minutes after the record of the fringe, i.e. correlated flux. 
As the atmospheric conditions change between both recordings, 
an additional error bar on the calibrated visibilities must be 
considered. Taking a similar error into account for both the 
SCI-PHOT and HIGH-SENS modes implicitly assumes that 
the variation in the transfer functions is much greater than 
the intrinsic error of one measurement due to the changing 
weather conditions. During good atmospheric conditions, this 
assumption is reasonable. During bad atmospheric conditions 
(i.e. up to a certain seeing threshold), the large fluctuation 
of the positions of the two beams reduces the beam overlap 
leading to strong variations in the fringe amplitude and the 
photometry. A solution to quantify the influence of the seeing on 
an individual measurement is to use the SCI-PHOT mode data 
with additional photometry (pseudo HIGH-SENS mode) taken 
under similar atmospheric conditions. Comparing calibrated 
visibilities reduced with the pure SCI-PHOT mode with those 
reduced with the pseudo HIGH-SENS mode gives a direct 
handle on the midterm (i.e. time span from fringe tracking to the 
photometry) influence of the seeing. Under bad conditions (i.e. 
seeing>1.5 "), the pseudo HIGH-SENS mode gives systematic 
lower visibilities, up to a level of 4%, when comparing them 
to the ones deduced from the pure SCI-PHOT mode, while the 
difference falls to zero when the seeing is smaller than 1.1". 
This is reasonable considering that the overlap of two beams is 
statistically reduced leading to systematically lower visibilities. 
We therefore consider an additional positive deviation of 4% 
(i.e. up to the 5% SCI-PHOT positive error) on the error bars 
associated to the calibrated visibilities recorded in HIGH-SENS 
mode under bad atmospheric conditions (i.e. seeing>1.5 "). 
Finally, the error bar on the calibrated visibilities recorded in 
HIGH-SENS mode under better atmospheric conditions (i.e. 
seeing<l.l ") is similar to the standard 10% error bar assigned 
to the SCI-PHOT mode measurements. 

The data quality has been checked using the MIA Graphical 
User Interface, and the software developed by Paris Observatory 
for MIDI 3 . Among the MIDI measurements, data sets #2, #5, 
#8, #15, and #25 (see Tables 2 and 3) led to outliers. Data 
sets #2, #15, and #25 show calibrated visibilities greater than 
unity, data set #5 shows an unusual spectral signature of the 
dispersed visibility, whereas data set #8 shows an abnormally 
high level of visibility amplitude compared to the other visi- 
bility data at close base, PA, and phase. Finally data sets #9, 
#10, #11, #12, #13, #14, and #16 show very broad dispersed 
power spectral densities (PSD) and/or fringe histograms (FH), 
especially from 8 to 10 pm since shorter wavelengths are more 
affected by the atmospheric conditions. If the PSD and FH are 
very broad, a very significant fraction of the fringe power is 
distributed outside the integration range, leading to a systematic 
underestimation of visibility. All the corresponding visibility 
measurements were discarded. 



2 http : //www . mpia-hd . mpg . de/MIDISOFT/ , 
http : / /www . strw . leidenuniv . nl/~nevec/MIDI/ 



3 Software is available through the JMMC website: 
http : //mariotti . u j f -grenoble . f r/ 
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Fig. 1. K-band uv coverage of the VINCI observations (diamonds) and 
N-band spectrally-dispersed (from 8 to 13 /jm) uv coverage of the MIDI 
observations (plus) (see Tables 1, 2, and 3). 
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Fig. 2. Broadband VINCI visibilities of R Scl (black diamonds) as a 
function of the projected baselines observed over 2 nights between 
November and December 2001 (see Table 1). Uniform disk diameters 
(blue squares to be read from the scale on the right axis) are computed 
from each visibility measurement. 



Figure 1 shows the K-band uv coverage of the VINCI 
observations (see Table 1) and the N-band spectrally-dispersed 
uv coverage of the MIDI observations (see Tables 2 and 3). The 
spatial coverage spreads over a large scale (from 15 to 60 m 
baselines), whereas the angular coverage is limited to 90°. 

Figure 2 presents the VINCI calibrated visibilities observed 
over 2 nights between November and December 2001 (see 
Table 1). The angular diameters of the equivalent uniform 
disk are computed from each visibility measurement and show 
a stability from phase-to-phase at different position angles. 
Figure 3 shows the MIDI calibrated visibilities observed over 14 
nights between January 2005 and December 2006 (see Tables 2 
and 3). The diameters of the equivalent uniform disk were 
computed from the visibilities at each spectral channel. The 
dimension of those diameters regularly increases from 9 to 13 
//m, indicating an extended circumstellar environment (Ohnaka 
et al., 2005). 

The phase shifts or differential phases correspond to the in- 
formation extracted from the difference between the phase at a 
given wavelength and the mean phase determined in the full N- 
band region. It provides information about possible asymmetry 
of the source in the region covered by the different interfero- 
metric measurements. As an example, strong phase signals were 
detected in sources harboring dusty disks (Deroo et al., 2007; 
Ohnaka et al., 2008). Due to the low values (< ±5°) of the cali- 
brated differential phases for each projected baseline, the global 
centrosymmetry of the close circumstellar environment of the 
object is confirmed in the mid-infrared for the corresponding an- 
gular coverage. 

3. Observed variability of RScl 

Sections 3.1, 3.2, and 3.3 discuss the variability of the star de- 
rived from photometry, polarimetry, and spectroscopy. Finally, 
Sect. 3.4 is dedicated to the interferometric variability giving 
constraints on the dynamic process at work in the close circum- 
stellar environment of R Scl. 



3.1. Photometric variability 

The K-band light curve of R Scl was discussed by Whitelock 
et al. (1997). The authors distinguish two main periods of 
374 and 1804 days. The shorter period is identified to be the 
pulsation period of the star, whereas the longer one might be 
a beat period between two similar pulsation periods (Houk, 
1963). More recent polarimetric measurements (Yudin & Evans 
2002, see Sect. 3.2) reveal a clumpy medium surrounding RScl 
where orbiting dusty clump structures could also be the origin 
of the long period. 

In the following, the phase of the star is established from the 
light curve according to 

(r-r )modP 

p ' C 1 ) 

where t corresponds to the observing time in Julian day, 
7o=2451 044 is the Julian date of the selected phase-zero point 
corresponding to the maximum light of the star (<f>o=0; see 
Fig. 4), and P=374 days is the pulsation period of the star. 

Besides the VINCI interferometric and MIDI spectro- 
interferometric observations of RScl, polarimetric (Yudin & 
Evans, 2002) and additional ISO/SWS spectrometric (Sloan 
et al., 2003) data of the star are available. 
Figure 4 represents the visual lightcurve of R Scl where the po- 
larimetric (Pola.), spectroscopic (ISO and MIDI), and interfero- 
metric (VINCI and MIDI) observing periods are marked. 

3.2. Polarimetric variability 

Wherever there is appreciable asymmetry in an astronomical 
situation, emitted light is likely to be polarized at some level. 
Yudin & Evans (2002) show that RScl displays polarimetric 
variability in the VRcIc bands on different time scales from 
hours to one year. The authors argue that the polarimetric 
variations of the star are most probably not related to binarity, 
since these variations occur on timescales that are inferior to one 
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Fig. 3. Spectrally-dispersed MIDI visibilities of R Scl (black thin lines to be read from the scale on the left axis) observed over 14 nights between 
January 2005 and December 2006 (see Tables 2 and 3). Uniform disk diameters (blue thick lines to be read from the scale on the right axis) are 
computed from the visibilities at each spectral channel. Error bars are not included for clarity. The data set numbers are defined in Tables 2 and 3. 
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Fig. 4. AAVSO visual lightcurve of RScl (Henden, 2009). Colored 
zones indicate the observing period of Polarimetric measurements (or- 
ange) (Yudin & Evans, 2002), ISO/SWS measurements (blue) (Sloan 
et al., 2003), VINCI measurements (violet), MIDI UT measurement 
(black), and MIDI AT measurements (green). The point surrounded by 
the open rectangle represents the chosen phase-zero point (0o=O). 



pulsation period (~374 days). With a wavelength-dependence of 
polarization close to Poc/l 4 , consistent with Rayleigh scattering 
by small («0.01 //m) amorphous grains, the authors suggest 



that the polarimetric variation would be related to dusty clumps 
located in the inner zone of the warm dusty shell (i.e. from 2 to 
4 R+) with a quasi-period of rotation in the range 470 to 1600 
days. 

The constraints provided by the MIDI differential phases 
(see Sect. 2), which do not exhibit any asymmetries in the 
close vicinity of the source (from 2 to 10 /?*), contradict the 
conclusions of Yudin & Evans (2002). Therefore, we are left 
with three possibilities. 

(i) The MIDI interfererometric angular coverage is 
outside the dusty clumps region. 

It is difficult to determine the position of the clumps because of 
the large uncertainty on their rotational period found by Yudin 
& Evans (2002). As it is not possible to locate the clumps at 
the time of the MIDI observations, we cannot exclude that the 
clumpy region is situated out of the interferometric coverage, 
which only spreads over an angle of 90° (see Fig. 1). 

(ii) The dynamic range of MIDI is not high enough to 
separate the emission of the dusty clumps from the one 
coming from the circumstellar material. 

We can investigate this point by an estimation of the constraints 
provided by the MIDI differential phase measurements. For 
that, we performed a perturbation of the intensity distribution 
of the best-fitting model (see Sect. 4.2) by a Gaussian-like dust 
cloud of 1.8 mas width (defined from the density structure found 
by Yudin & Evans 2002) varying its flux contrast, its distance 
(from 2 to 4 and its position angle. Our tests allow any 
mid-IR dusty clump source representing more than about 4% 
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of the total flux of the object (i.e. 4 to 8 Jy between 8 pm and 
13 pm) to be discarded. Taking the condensation density p c and 
the condensation mass derived by Yudin & Evans (2002) into 
account, and describing the dusty clump as a blackbody of 1000 
to 1500 K at a distance ranging from 2 to 4 R+, we obtain a flux 
range between 2.0 and 3.5 Jy in the N-band. As these values 
are below 4% of the total flux of the object, this solution can be 
retained. 

(iii) The dusty clumps emission vanished from the N- 
band spectral window during the 10 years separation from 
the polarimetric to the MIDI interferometric observations. 

Shock waves could push the dusty clumps observed by polarime- 
try to regions where its emission peak is outside the mid-infrared 
window. However, taking a shell expansion velocity from 10.5 to 
16.5 kms" 1 into account for RScl (Schoier & Olofsson, 2001; 
Wong et al., 2004), the dusty clump would be located around 
14.4 to 23.4/?* (within the FOV of both AT and UT telescopes) 
at the end of the 10-year separation from the polarimetric to the 
MIDI observations. Assuming a Planckian radiation field geo- 
metrically diluted with distance from the star, the temperature 
of the clump reaches about 400 K corresponding to an emission 
peak well inside the mid-infrared window. Therefore this solu- 
tion can be discarded. 

3.3. Spectrometry variability 

As discussed in Sect. 2, most of the mid-infrared flux originates 
in the inner 300 mas of the source. Therefore, the mid-infrared 
spectra obtained from the AT and UT single dishes are di- 
rectly comparable to the ones produced by the ISO telescope 
(FOV=33"x20")- To keep consistency in the comparison to the 
measurements, only ISO/SWS spectra using the AOT 1 template 
(Leech et al., 2003) are considered in the following. Figure 5 
compares both calibrated MIDI and ISO/SWS spectra (Sloan 
et al., 2003) of R Scl around the maximum visual light (0~O). 
Figure 6 presents the calibrated MIDI and ISO/SWS flux of 
RScl as a function of the phase at 8.5 pm (pseudo-continuum: 
defined where the blackbody continuum fits to the ISO/SWS 
spectrum of the star), 11.3 pm (SiC feature), and 12.5 /mi 
(C 2 H 2 +HCN) spectral bands (Hron et al., 1998; Gautschy-Loidl 
et al., 2004; Zijlstra et al., 2006). Each monochromatic flux is 
determined by integration of the narrow-band filter fluxes over 
a bandwidth of AA-0.2 pm. The MIDI error bars are estimated 
from the standard deviation of the N-band flux-calibrated 
spectra determined from each night and each telescope, i.e. 
typically 3 to 10 spectra, except at phases 0.23 and 0.79 where 
only one spectrum is available (see Tables 2 and 3). 

From Figs. 5 and 6, we can argue that the level of flux re- 
ceived by MIDI is equivalent to the one of ISO. It means that the 
cycle-dependent flux level of the object in the N-band has prob- 
ably not changed significantly from the ISO/SWS (1996/1997) 
to the MIDI (2005/2006) observations. From Fig. 6, it is rather 
difficult to determine a real trend in variability, keeping in mind 
typical errors in the absolute ISO/SWS flux level ranging from 
10 to 12% (van der Bliek et al., 1998) and an error of about 20% 
for the MIDI flux. 

3.4. Interferometric variability 

To study the interferometric variability of the star, we compare 
visibility measurements observed at the same projected baseline 
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Fig. 5. Comparison of the ISO/SWS spectrum of R Scl at phase 0.97 
(thick line), and the MIDI spectrum of R Scl at phase 0.98 (error bars). 




Fig. 6. MIDI (squares) and ISO/SWS (triangles) flux of R Scl as a func- 
tion of the visual phase in the 8.5 (left), 11.3 (middle), and 12.5 /im 
(right) spectral band corresponding to the pseudo-continuum, the SiC 
feature, and the C2H2+HCN feature, respectively. 



lengths for very close position angles (few degree differences) 
from one epoch to another one. The same projected baselines 
need to be used to probe the same spatial regions of the object, 
whereas very close projected angles are needed to avoid effect 
coming from potential deviation in the geometry of the object 
from spherical symmetry. 

We first compare both sets of VINCI visibility measure- 
ments at 0=0.17 and 0=0.23 for the given 16 m projected 
baseline length through a range of 5° in position angle. This 
comparison does not display any significant (larger than 2cr) 
temporal variation. This fact is, however, not surprising given 
the short time interval at post-maximum brightness of the star 
and because the broadband measurements will average out the 
possible short term effects of shock fronts in CO or CN lines 
(e.g. Nowotny et al. 2005a, Paladini et al. 2009). As the VINCI 
data only probe the very low spatial frequency, what remains 
are thus only the slow variations due to stellar pulsation (P=374 
days). 

In the case of the MIDI measurements, the lowest spectral 
resolution (R=30) is high enough to resolve the broad features 
contained in the N-band spectral window (i.e. C2H-2, HCN, C3, 
and SiC). Therefore, we do not expect any effect that could pos- 
sibly smear out the temporal variation coming from emission or 
absorption bands in that region. Figure 7 shows comparisons be- 
tween the MIDI spectrally-dispersed visibilities from phase-to- 
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Fig. 7. Comparison between spectrally-dispersed MIDI visibilities from 
phase-to-phase (black thin dashed and dotted lines) at a given projected 
baseline length for very close position angles. The data set numbers are 
defined in Table 3. 

phase at a given projected baseline length for very close position 
angles. The measurements do not show any significant variabil- 
ity effect for both 15 and 31 m baselines. Because those baselines 
probe the region located from ~17. 5 to 100 mas (3.4 to 19.6R+), 
this means that the stellar radiative pressure is not strong enough 
to reveal a significant movement of the warm mass dust shells 
detectable at the spatial resolution of the interferometer (from 
~8 to 12 mas: 2.9 to 4.3 AU at 360pc). Such behavior agrees 
with the amplitude of the mass shell of theoretical carbon-rich 
dynamic models, which does not exceed 3 AU over a full cycle 
in the mid-infrared (see Paladini et al. 2009). 

4. Modeling the dynamic atmosphere of R Scl 

The strategy of determining the best-fitting parameters for 
the dynamic model atmosphere of the C-rich star RScl is the 
following: 

Step 1: The hydrodynamic computation (see Sect. 4.3.1) 
starts with an initial hydrostatic structure mainly characterized 
by the following parameters (Hofner et al., 2003): luminos- 
ity, mass, effective temperature, metallicity, C/O ratio, and 
micro-turbulent velocity (for line opacities). In order to derive 
these parameters, we first perform a hydrostatic modeling of 



the stellar atmosphere (Sect. 4.1) using a grid of COMARCS 
models (Aringer et al., 2009). By fitting synthetic spectra and 
visibilities based on such models to observed near-infrared 
spectro-interferometric measurements of R Scl, we find the best 
set of stellar atmospheric parameters for the initial hydrostatic 
structure of the dynamic model. 

Step 2: The mid-infrared excess of RScl shows an SiC 
feature around 11.2 yum (see Fig. 5). This component must 
therefore be considered in the spectral synthesis for the com- 
parison with observational data. Since only AmC dust is treated 
in the hydrodynamic calculation (Hofner et al., 2003), opacities 
of SiC dust are added in the a posteriori radiative transfer 
computation (see Sect. 4.3.1). Those opacities are determined 
by finding the proper SiC optical constants, and the fractional 
abundance of SiC to AmC dust that has condensed. These 
two characteristics are determined in Sect. 4.2.1 by comparing 
observed mid-infrared spectro-interferometric measurements of 
RScl with synthetic spectra and visibilities derived from the 
DUSTY numerical code (Ivezic & Elitzur, 1997; Ivezic et al., 
1999), which solves the radiative transfer equation through a 
dusty envelope. 

In Sect. 4.2.2 we present the disagreement of the best- 
fitting model created by the combination of the hydrostatic 
COMARCS atmosphere and DUSTY circumstellar envelope 
(COMARCS +DUSTY model) with the observational data. This 
shows the necessity of using a self-consistent hydrodynamic 
model to understand how the time-dependent processes influ- 
ence the atmospheric structures of the star. 

Step 3: Finally, the self-consistent dynamic model used in 
Sect. 4.3.3 is chosen by the parameters deduced from the best 
hydrostatic fit of the stellar atmosphere (see Sect. 4.1) together 
with the SiC dust parameters derived from the DUSTY modeling 
(see Sect. 4.2.1). We then interpret the spectro-interferometric 
data of RScl using this model to discuss the dynamic picture 
deduced from the observations. 

4. 1 . Determination of the stellar parameters of R Scl with 
hydrostatic model atmospheres 

Besides the ISO/SWS spectra (see Sect. 3. 3), ground-based pho- 
tometric measurements (Bagnulo et al., 1998; Whitelock et al., 
2006) are available for R Scl (see Table 4). Using the relation 
of Milne & Aller (1980), also used by Bagnulo et al. (1998), a 
visual interstellar extinction value of 0.18 was derived. The ra- 
tio of total-to-selective extinction at V, 3.1, on average, allowed 
us to correct the other magnitudes using the absorption curve of 
Fitzpatrick (1999). The data of Whitelock et al. (2006) are dis- 
persed over a wide range of 26.3 pulsation periods of the star, 
whereas the data of Bagnulo et al. (1998) are measured at post- 
maximum light of the object. 

Because of all the absorption lines coming from strong 
molecular opacity sources (mainly species containing C, H, N, 
O; see Gautschy-Loidl et al. 2004), spectra of cool stars differ 
a lot from a blackbody continuum spectra as is highlighted in 
Fig. 8. The use of a more realistic synthetic carbon-rich spec- 
trum is then required. The COMA code, developed by Aringer 
(2000), solves for the ionization and chemical equilibrium 
(using the method of Tsuji 1973) at a given temperature and 
density (or pressure) combination for a set of atomic abundances 
assuming LTE. Based on these results, atomic and molecular 
opacities are calculated. Details of computation and information 
on recent updates on the COMA code can be found in Lederer 
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Table 4. Broadband photometry of R Scl 



Band 


Extinction 


Dereddened 


Phase 


Reference 




coefficient 


magnitude 






B 


0.25 


11.03+0.3 1 


0.17 


Bagnuloetal. (1998) 


V 


0.18 


4.55±0.7 1 


0.17 


_ 


J 


0.05 


1.56±0.05 


0.12 




H 


0.03 


0.46±0.05 


0.12 




K 


0.02 


-0.24±0.06 


0.12 




L' 


0.01 


-0.84±0.28 


0.12 




J 


0.05 


1.97±0.86 


Whole cycle 


Whitelock et al. (2006) 


H 


0.03 


0.63±0.64 






K 


0.02 


-0.10±0.36 






L 


0.01 


-0.74±0.28 







'Uncertainties calculated from the measurements given in the Simbad 
database. 



& Aringer (2009), and Aringer et al. (2009). These opacities 
are then introduced in a spherically symmetric radiative transfer 
code to calculate the emergent intensity distribution. This inten- 
sity distribution is then used to derive the synthetic spectrum 
and visibility profile, which are compared with observational 
data of the star. 

The determination of the stellar parameters is done by 
comparing hydrostatic COMARCS models to the spectropho- 
tometric data up to 5 yum. This limit is fixed because of the 
contribution of the dust becoming non negligible at longer 
wavelengths for optically thin, warm dusty objects like RScl 
(Kraemer et al., 2002). This will be confirmed by the DUSTY 
modeling in Sect. 4.2.1. The hydrostatic model was chosen as 
best-fitting of the broadband photometric data (see Table 4) 
and the various features in the ISO/SWS spectra (see Fig. 8) 
averaged over all the phases. Furthermore, the best-fitting model 
should also agree with the K-broadband VINCI interferometric 
measurements (see Sect. 2). 

A least-square fitting minimization was completed from the 
grid of models included in Aringer et al. (2009), varying the 
parameter values of effective temperature of the star from 2600 
to 3000 K (AT eff =100 K), the luminosity from 6000 to 8000 L 
(AL=200 L ), the distance from 300 to 500 pc (Ad=50pc), the 
mass for values of 1 and 2M , the metallicity for values of 
1, 0.33, and 0.1 Z©, and the C/O ratio for values of 1.1, 1.4, 
and 2.0. Finally, a value of 2.5 kms -1 for the microturbulent 
velocity (£), coming from the fit of very high spectral resolution 
data of carbon stars, is taken into account. 
We obtained the best-fitting model to the ISO/SWS spectromet- 
ric data from 2.4 to 5 //m with a reduced^ 2 of 3.5 (see Table 5). 

Figure 8 shows the best-fitting hydrostatic model on the 
spectrophotometric data of the star. As the VINCI data are 
integrated over the K-band (2.0-2.4 /mi) defined by the corre- 
sponding filter curve, the synthetic broadband visibility results 
in a superposition of all the monochromatic profiles over the 
width of the band (Kervella et al., 2003). Figure 11 presents 
the best-fitting hydrostatic model superimposed on the VINCI 
visibility measurements. 

Table 5 summarizes the stellar atmospheric parameters of 
RScl deduced from the best-fitting hydrostatic model of the 
spectrophotometric and near-infrared interferometric measure- 
ments. 
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Fig. 8. Best-fitting hydrostatic model (blue dashed line) on the spec- 
trophotometric measurements. Photometric measurements (squares: 
Bagnulo et al. 1998 + triangles: Whitelock et al. 2006; see Table 4) and 
ISO/SWS spectra of R Scl at different epochs (dotted lines) are shown. 
A blackbody continuum having the same temperature (T cff =2700 K) 
and angular diameter (0^ = 10.2 mas) as the synthetic COMARCS spec- 
trum is superimposed (black solid line). 

Table 5. Stellar atmospheric parameters of R Scl deduced from best- 
fitting hydrostatic model of the spectrophotometric and near-infrared 
interferometric measurements. 



Parameter 


Value 


Effective temperature: T e g (K) 


2700 


Luminosity: L (L Q ) 


7000 


Distance: D (pc) 


350 


Central star diameter: 0* (mas) 


10.2 


Surface gravity: log g 


-0.7 


Micro-turbulent velocity: f (kms~') 


2.5 


Stellar mass: M (M ) 


1 


Metallicity: Z (Z Q ) 


1 


Carbon over oxygen ratio: C/O 


1.4 



The values found agree well with the ones found in the liter- 
ature. The distance of 350 pc is relatively far from the Hipparcos 
distance of 474 pc but in very good agreement with the distance 
of 360 pc derived from the period-luminosity relation for car- 
bon stars (Groenewegen & Whitelock, 1996; Whitelock et al., 
2008). The angular diameter of the central star is relatively close 
to the one coming from the bolometric luminosity of R Scl (from 
11.8 to 13.8 mas; Yudin & Evans 2002). The C/O ratio of 1.34 
derived by Lambert et al. (1986) agrees with the value that we 
found. The surface gravity, mass, metallicity, and C/O ratio are 
also in good agreement with those used by Hron et al. (1998) 
(logg=-0.5/-0.6, M=l M , Z=l Z , and C/0=1.4, respectively). 

4.2. Including of a dusty environment 

4.2.1 . Determination of the SiC dust parameters 

The DUSTY code (Ivezic & Elitzur, 1997; Ivezic et al., 1999) 
solves the problem of radiation transport in a circumstellar dusty 
environment by calculating the radiative transfer equation in 
plane-parallel or spherical geometry. The radiation from the star 
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is scattered, absorbed, and re-emitted by dust grains. DUSTY 
has built-in optical properties for the most common types of 
astronomical dust. The stellar spectrum can be specified in a 
numerical form as a separate user-supplied input file. In our 
case, the best-fitting COMARCS stellar atmosphere model, 
defined in Sect. 4.1, is used as central source. The code can 
compute the wind structure by solving the hydrodynamic equa- 
tions, including dust drift and the stars gravitational attraction, 
as a set coupled to radiative transfer (Ivezic et al., 1999). The 
calculation is performed for a typical wind in which the final 
expansion velocity exceeds 5 km s _1 , in agreement with the 
expansion velocity found for RScl from 10.5 to 16.5 kms" 1 
(Schoier & Olofsson, 2001; Wong et al., 2004). This allows 
derivation of the mass-loss rate of the object, fixing the dust 
grain bulk density p s and the gas-to-dust mass ratio r g d. 
The input parameters of DUSTY related to the dusty circum- 
stellar environment are the optical constants and the fractional 
abundances of the relevant dust grains, the grain size distri- 
bution, the overall radial optical depth at a given wavelength, 
the geometrical thickness, and the dust temperature at the 
inner boundary. Different types of grains can have different 
temperatures at the same location. However, DUSTY currently 
treats mixtures as single-type grains whose properties average 
the actual mix, therefore only one temperature is specified. 

The determination of the parameters is done by fitting the 
COMARCS +DUSTY model to the photometric measurements, 
the ISO/SWS spectrometric data from 0.4 to 25 /mi, and the 
MIDI spectrometric measurements from 8 to 13 /mi, aver- 
aged over all the phases. Furthermore, the best-fitting model 
should also agree with the K-band VINCI and N-band MIDI 
interferometric measurements. For that, a least-square fitting 
minimization was completed first by testing different sets of 
optical constants for AmC and SiC grains and by computing a 
large grid of models varying the fractional abundances of the 
relevant dust grains ranging from 100% AmC to 100% SiC 
grains in 10% increments, the 0.55 //m overall radial optical 
depth from 0.05 to 1 (Aro.ss^m^O.OS), and the dust temperature 
at the inner boundary from 800 to 1500 K (Ar in =100 K). 
As the geometrical thickness does not have any significant 
impact on the model spectrum because of the very optically 
thin circumstellar environment (Tn.2^m=7.6xl0 -3 , see Table 6), 
we fixed this parameter to 1000 dust shell inner boundary radii 
knowing that the result is similar for values between 100 to 10 4 . 
The grain size distribution was also fixed to a standard MRN 
grain size distribution as described by Mathis et al. (1977) 4 . 
We obtained the best fit model to the ISO/SWS spectromet- 
ric data from 2.4 to 25 fim with a reduced^ 2 of 3.5 (see Table 6). 

Figure 9 shows the best-fitting COMARCS +DUSTY model 
on the spectrophotometric data of the star. Figures 10 and 11 
present the best-fitting COMARCS +DUSTY model superim- 
posed on the MIDI and VINCI visibility data, respectively. 
Table 6 summarizes the parameters of the dusty envelope found 
from the best-fitting COMARCS +DUSTY model on the spec- 
trophotometric and the near- and mid-infrared interferometric 
measurements. These parameters are complementary to the 
stellar atmospheric parameters (see Table 5). 

The resulting values agree well with the ones found in the 
literature. The inner radius of the dust shell is close to 5 R+ (see 



4 n(d) oc a 3,5 with minimum and maximum grain sizes of 0.005 jim 
and 0.25 fim, respectively 
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Fig. 9. Best-fitting COMARCS+DUSTY model (red solid line) on the 
ISO/SWS (dotted lines), MIDI (error bars) spectrometric measure- 
ments, and photometric (squares and triangles; see Table 4) measure- 
ments of R Scl. The stellar contribution (blue dashed line) and the dust 
shell contribution (green dashed-dotted line) are also shown. 
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Fig. 11. Best-fitting hydrostatic model (blue dashed line) and best- 
fitting COMARCS+DUSTY model (red solid line) on the VINCI visi- 
bilities of R Scl (error bars). 



Table 5) as determined by Lorenz-Martins et al. (2001). The op- 
tical constants of AmC from Rouleau & Martin (1991) and SiC 
from Pegourie (1988), together with the fractional abundance 
of SiC/AmC of 0.08 found by Lorenz-Martins et al. (2001), 
are also in good agreement with the values found in this work. 
These last two characteristics allow us to determine the opac- 
ities of SiC used to interpret the spectro-interferometric data 
of RScl with the help of the dynamic model (see Sect. 4.3.3). 
The AmC dust grain optical constants from Rouleau & Martin 
(1991) are the same as adopted for the hydrodynamic calcula- 
tion (Andersen et al., 2003; Hofner et al., 2003). Good agree- 
ment is also found with the works of Dehaes et al. (2007) and 




Fig. 10. Best-fitting COMARCS+DUSTY model (solid line) on the spectrally-dispersed MIDI visibilities of RScl (error bars) from the smallest 
to the largest projected baselines. The data set numbers are defined in Tables 2 and 3. 



Table 6. Dusty circumstellar envelope parameters of RScl deduced 
from the best-fitting COMARCS+DUSTY model on the spectropho- 
tometric and the near- and mid-infrared interferometric measurements. 



Parameter Value 



Shell inner radius: e in (mas) 


22.4 


Inner boundary temperature: Tj n (K) 


1200 


Grain chemical composition 1 


90% AmC + 10% SiC 


Grain size distribution 2 


MRN 


Geometrical thickness: Y (inner radii) 


1000 


Visual optical depth: T 5 5llm 


0.40 


1.0 /im optical depth: Ti. 0/Jm 


0.18 


11.2 yum optical depth: T ll 2llm 


7.6xl0~ 3 


Mass-loss rate 3 : M (M yr _1 ) 


(6.7+2)xl0- 7 



'AmC stands for amorphous carbon of Rouleau & Martin (1991) and 
SiC for silicon carbide of Pegourie (1988). 
2 grain size distribution as described by Mathis et al. (1977). 
'assuming a dust grain bulk density p s = 1 .85 g cirr 3 (Rouleau & Martin, 
1991) and a gas-to-dust mass ratio r gd =590 (Schoier et al., 2005). 



Lorenz-Martins et al. (2001), who determined a visual opti- 
cal depth of 0.5^ 3 and a 1.0 /im optical depth of 0.1, respec- 
tively. The dust shell inner boundary temperature of lOOO^oo K 
found by Dehaes et al. (2007) agrees well with the value deter- 
mined in this work. Finally, the mass-loss rate of the star that we 
found also well agrees with the ones determined by Le Bertre 
(1997): M=10~ 7 M yr'; Gustafsson et al. (1997): M=4xl0~ 7 
M yr" 1 ; and Wong et al. (2004): M=2 to 5xl0~ 7 M Q yr" 1 . Such 
low values of mass-loss rate are consistent with the very opti- 
cally thin circumstellar environment that we found from the fit 

(Tll. 2/ ,m=7.6xl0- 3 ). 



4.2.2. Comparison of the COMARCS+DUSTY model with 
observations 

While the best-fitting COMARCS+DUSTY model is able to 
reproduce the largest part of the ISO/SWS spectra (see Fig. 9), 
we can see a discrepancy from 5 to 7 /mi, which could be 
linked to wrong or incomplete C3 line list opacities in that 
spectral range 5 . Another discrepancy between the model and 
the ISO/SWS spectro metric data appears around the 13.7 //m 
C2H2 absorption feature. Gautschy-Loidl et al. (2004) argue 
that this weak and sharp 13.7 /im absorption would come 
from cool layers above the 500 K dust and gas layer proposed 
by J0rgensen et al. (2000) to explain the lack of the 14 //m 
photospheric band. 

The model matches the observed near- and mid-infrared 
visibilities data rather well (see Figs. 10 and 11), except for 
few measurements in the mid-infrared. The first disagreement 
comes from the amplitudes of the model visibilities, which 
are either higher or lower than the amplitudes of the visibility 
measurements at certain baselines. The most probable reasons 
for that follow. 

(i) The geometry of the object: 

Potential deviations from spherical symmetry (corresponding to 
the geometry of the COMARCS+DUSTY model) of the object 
at certain position angles could explain the disagreement. As 
the object subsequently appears more or less extended from 
one position angle to another, the corresponding amplitude 

5 Since there is only one C3 line list available, a cross check of dif- 
ferent line lists is not possible. 
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of the visibility measurements decreases or increases, respec- 
tively, in contrast to the constant visibility amplitude of the 
COMARCS+DUSTY model. 

(ii) The variability of the object: 

It has been shown that the motion of the dynamic model atmo- 
spheres without mass loss mainly follows the pulsation of the 
stellar interior, whereas in the case of mass-losing models, the 
dust formation may deform the corresponding sinusoidal pattern 
(see Fig. 2 of Hofner et al. 2003). As the object subsequently 
appears more or less extended from phase-to-phase, the corre- 
sponding amplitude of the visibility measurements decreases 
or increases (Paladini et al., 2009), respectively, in contrast to 
the constant visibility amplitude of the COMARCS+DUSTY 
model. 

Another phenomenon that cannot be reproduced by the 
model is the slope in the N-band visibility data for the largest 
60 m projected baselines in the 8-9 pm spectral range (see 
Fig. 10). As the model overestimates the visibility measure- 
ments, this reveals that the model is too optically thin. The 
regions probed at these baseline and spectral band are located 
between the photosphere of the star (~5 mas; see Table 5) and 
the dust shell inner radius (~22.4mas; see Table 6). Because 
of the very low optical depth of the model in the 8-9 /im 
spectral range [f(8-9 yum)=4.6xl0~ 3 ], leading to domination 
of the stellar flux over the flux of the dust shell, the model is 
therefore close to the central star's intensity distribution. A way 
to reproduce the MIDI measurements is to increase the optical 
thickness of the model at the corresponding spatial frequencies. 
This can be done by adding extended molecular layers of C2H2 
and HCN emitting in that spectral range and located above the 
stellar photosphere. 

Many works in interferometry revealed that diameters of Mira 
and non-Mira M- and S-type stars appear systematically larger 
than expected in the near-infrared (Mennesson et al., 2002; 
Ohnaka, 2004; Perrin et al., 2004; Weiner, 2004; Verhoelst et 
al., 2006) and in the mid-infrared (Weiner et al., 2000, 2003a,b; 
Weiner, 2004; Ohnaka, 2004; Ohnaka et al., 2005; Verhoelst et 
al., 2006; Sacuto et al., 2008). This increase cannot be attributed 
only to dust shell features, but also to the possible existence of 
extended gaseous shells above the stellar photosphere. Such a 
molecular shell, also called MOLsphere (Tsuji et al., 2000), is 
favored by the low surface gravity (log g--0.7) and the high 
luminosity (~7000 L Q ) of the star, allowing the levitation of 
the upper layers by shock waves or other dynamic phenomena 
in an easiest way. The temperatures in this levitated matter are 
low enough that large amounts of polyatomic molecules can 
form (Hron et al., 1998; Woitke et al., 1999; Gautschy-Loidl et 
al., 2004; Verhoelst et al., 2006). Ohnaka et al. (2007) present 
their results on the modeling of the MIDI measurements of 
the carbon-rich Mira variable VOph. The model consists of 
optically thick warm and cool C2H2 molecular layers and an 
optically thin AmC+SiC dust shell. 

Although this MOLsphere is useful for increasing the 
optical thickness of the close circumstellar structure, it remains 
an inconsistent ad-hoc model. In the classical MOLsphere 
scenario, as is usually applied to interpreting interferometric 
or spectroscopic data, it is assumed that there are one or a few 
static layers with a certain concentration and temperature of a 
molecular species. The possible physical status and origin of 
such a layer is left unanswered, which is problematic for the 
following points. 



(i) The radial temperature-density distribution of gas and dust 
around evolved objects is always a continuous and inconstant 
phenomenon. 

(ii) The molecular (and dust) concentrations in a layer depend 
on each other, so it is impossible to freely choose the amount of 
certain species. 

(iii) The radial structures are extremely time dependent. 

As shown by Hron et al. (1998), Woitke et al. (1999), and 
Gautschy-Loidl et al. (2004), the hydrodynamic modeling leads 
to a self-consistent formation of extended molecular structures 
located between the photosphere and the dust formation zone. 
Finally, Hofner et al. (2003) argue that sequences of hydrostatic 
models with varying stellar parameters cannot reproduce the ef- 
fects of shock waves, levitation, or winds. The only way to un- 
derstand how the dynamic processes influence the atmospheric 
structure on different spatial scales is the use of time-dependent 
dynamic atmospheric modeling. 

4.3. Dynamic model atmosphere for RScl 

The aim of this section is to show that a self-consistent hydro- 
dynamic modeling forming a gaseous and dusty environment by 
itself is able to reproduce the time-dependent signature of the 
spectro-interferometric data of a carbon-rich, long-period vari- 
able star. Section 4.3.1 presents the dynamic model atmosphere 
used for the interpretation. Section 4.3.2 connects the bolometric 
phases given by the model with the visual phases of the observa- 
tions which is needed to compare time-dependent models with 
real measurements of RScl. Section 4.3.3 presents the specific 
dynamic atmospheric model that we chose for the comparison 
with the spectro-interferometric data of the star. Section 4.3.4 
shows the ability of the dynamic model to form molecular lay- 
ers located above the photosphere of the star in a self-consistent 
way. Finally, the last section (Sect. 4.3.5) compares the dynamic 
model with the time-dependent spectro-interferometric data. 

4.3.1. Presentation of the hydrodynamic model atmospheres 

The self-consistent hydrodynamic structures for dust-forming 
atmospheres and winds used here are taken from the models 
contained in Mattsson et al. (2008) and from additional mod- 
els required for the present work (Mattsson, priv. comm.). The 
dynamic model atmosphere is based on the code described by 
Hofner et al. (2003) and Mattsson et al. (2010). The models were 
obtained by solving the coupled system of hydrodynamics and 
frequency-dependent radiative transfer in spherical symmetry, 
together with a set of equations describing the time-dependent 
formation and the growth and evaporation of dust grains. The 
formation of AmC dust is treated in a self-consistent way us- 
ing the moment method described by Gail & Sedlmayr (1988) 
and Gauger et al (1990). Optical properties of AmC grains were 
taken from Rouleau & Martin (1991). Each dynamic model 
starts from a hydrostatic initial structure, and the effects of stel- 
lar pulsation are simulated by a variable inner boundary just be- 
low the stellar photosphere (piston model). The dynamic code 
provides temporally varying temperature and density stratifica- 
tions of the object. Based on such radial structures, at given in- 
stances of time, opacities are computed using the COMA code 
(see Sect. 4.1). The resulting opacities are used as input for a 
spherical radiative transfer code that computes the emergent in- 
tensity distribution for every frequency point of the calculation. 
This is used to calculate the synthetic spectra and visibility pro- 
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files, which are compared to the observational data of RScl. 
Since only AmC dust is treated in the hydrodynamic calcula- 
tion, opacities of SiC dust are added in the a posteriori radia- 
tive transfer computation. Those opacities are determined with 
COMA taking SiC optical constants into account and assuming 
that SiC dust is scaled with AmC dust. The scaling factor corre- 
sponds to the fractional abundance of SiC to AmC dust, which 
has condensed. Finally, SiC dust is assumed to follow the same 
temperature-density structure as AmC dust, derived from the full 
dynamic calculation. 

4.3.2. Bolometric and visual phases 

The phases of the dynamic model correspond to the bolometric 
phases derived from its bolometric lightcurve (Nowotny et al., 
2005b). The phase-zero point was chosen at the maximum bolo- 
metric luminosity of the given dynamic model. As the radiation- 
hydrodynamic code stores atmospheric structures at more or less 
randomly distributed instances of time, the uncertainty in the 
choice of the phase-zero point amounts to a range of +0.03. 
For O-rich Miras, a small shift between bolometric and visual 
lightcurves probably exists, such that bolometric phases lag be- 
hind visual ones by a value of w 0.1 (Lockwood & Wing, 1971; 
Nowotny et al., 2010). This phase shift is caused by the behav- 
ior of molecules occurring in atmospheres of red giants and their 
corresponding spectral features in the visual. Recently, Nowotny 
et al. (2010) have studied the photometric variation in the V-filter 
for similar dynamic model atmospheres as the ones used here. 
The authors find that, for models developing no wind such as 
their model W or the one used in this paper to represent R Scl 
(see Sect. 4.3.3), the visual lightcurves follow the bolometric 
light variation caused by the variable inner boundary (piston) 
closely. Therefore, in the following, bolometric phases of the 
model within a range of ±0.03 around the observed visual phases 
of R Scl are used for comparisons between model and observa- 
tions. 



4.3.3. Choice of a specific dynamic model 

The choice of the best dynamic model for R Scl comes from the 
stellar parameters and the distance derived from the previous 
hydrostatic modeling (see Sect. 4.1). 

Among the dynamic models that we have, one model shows 
parameter values of the hydrostatic initial structure close to the 
ones we derived from the hydrostatic COMARCS modeling 
(see Table 5). Table 7 compares the stellar parameters of the 
hydrostatic model initiating the hydrodynamic computation with 
the stellar parameters derived from the best-fitting hydrostatic 
model determined in Sect. 4. 1 . 

Opacities of SiC dust are added in the a posteriori radiative 
transfer calculation using the SiC optical constants of Pegourie 
(1988) and a fractional abundance of SiC to AmC dust of 
10%, as derived from the COMARCS +DUSTY modeling (see 
Sect. 4.2.1). 

Two other free parameters are required for the hydrody- 
namic computation (see Sect. 4.3.1): 

P mo d: the pulsation period of the piston. 

A pulsation period of P mo d=390 days, the closest to the pulsation 



Table 7. Parameter values of the hydrostatic initial structure of the hy- 
drodynamic model used in this work (first row) compared to the pa- 
rameters derived from the best-fitting hydrostatic model determined in 
Sect. 4.1 (second row). 



Model 


U 




T* 


z* 


<r 


C/O 




(Lq) 


(Mo) 


(K) 


(Z ) 


(kms- 1 ) 




Hydrodynamic 


7080 


1 


2800 


1 


2.5 


1.35 


Hydrostatic 


7000 


1 


2700 


1 


2.5 


1.40 



period of R Scl (374 days), was used. 

Au p : the amplitude of the piston velocity. 

Different amplitudes of the piston velocity were tested from a 
range of Au p =2 to 6 kms -1 , with steps of 1 kms -1 (except for 
the Au p =3 km s _1 model not at our disposal). 
Determination of the best-fitting piston velocity amplitude is 
done by comparing the synthetic spectra based on these models 
to the spectrophotometric data. Increasing the amplitude of 
the piston moves the material to larger distances from the 
star, which favors dust formation and reduces the effective 
gravitation (e.g. Winters et al. 2000). 

For a piston velocity lower than 4 km s _1 (i.e. model with 
Au p =2 kms -1 in Fig. 12), it appears that the model shows 
too low a level of flux compared to the spectrometric data. 
Furthermore, as already discussed in Gautschy-Loidl et al. 
(2004), a broad and intense C2H2/HCN absorption band around 
14 jum is present, in pronounced contrast to the observed 
ISO/SWS spectra which show a continuum-like infrared excess 
with a weak and narrow C2H2 feature at 13.7 /im (see Fig. 12). 
Finally, the variability amplitude of the model is smaller than 
0.3 in the V magnitude, as compared with the 1.4 variability 
amplitude of the cycle-to-cycle averaged visual magnitudes of 
the star. Therefore, in the case of too low a piston velocity, the 
dynamic model corresponds to a dust-free pulsating atmosphere 
without any wind (no mass-loss), leading to a pure photospheric 
spectrum such as the one shown by the blue dashed line of 
Fig. 9. 

When the piston velocity amplitude exceeds 4 km s 1 (i.e. 
model with Au p =6 km s _1 in Fig. 12), the dynamic model shows 
too high a level of flux compared to the spectrometric data and 
a filling up of the whole molecular features at wavelengths 
longwards of 3.5 jum. Furthermore, the variability amplitude of 
the model is greater than 5.5 in the V magnitude, as compared 
with the 1 .4 variability amplitude of the cycle-to-cycle averaged 
visual magnitudes of the star. In that case, the dynamic model 
develops a strong wind (M>10~ 6 M yr _1 ) where the large 
amount of formed gas and dust, located in the outer layers, 
obscures the innermost structures of the stellar atmosphere, 
although not significantly enough to fill up the 3.05 fim feature. 
A dynamic model with a piston velocity of 4 km s _1 gives the 
best overall match of the spectrophotometric data from 0.4 to 
25 /im, even if the emission is still lower than the [13-25 /im] 
infrared excess of the ISO/SWS spectrum (see Fig. 12). The 
variability amplitude of the model of about 1 in the V magni- 
tude is relatively close to the 1.4 variability amplitude of the 
cycle-to-cycle averaged visual magnitudes of the star. In that 
case, the model is, in some sense, a transition case between 
a simple, periodically pulsating dust-free atmosphere and an 
object featuring a dusty wind. The kinetic energy supplied by 
the pulsation (piston) causes strong, quasi-periodic atmospheric 
dynamics and some intermittent dust formation, but is not 
sufficient to start a dust-driven outflow. 
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Fig. 12. Synthetic spectra based on the dynamic models at phase 0.98 
with piston velocities of 2 (green line), 4 (red line), and 6 km s _I 
(blue line), keeping the same parameter values as the ones contained 
in Table 7. The dotted line is the ISO/SWS spectrum of the star at 
phase 0.97, whereas the error bars in the N-band correspond to the MIDI 
fluxes at phase 0.98. Squares and triangles correspond to the photomet- 
ric measurements of the star (see Table 4). 



The possible reasons for the discrepancy in the 5 to 7 //m wave- 
lengths regions are similar to the ones expressed in Sect. 4.2.2 
suggesting wrong or incomplete C3 line list opacities in that 
spectral range. The model is also unable to reproduce the weak 
observed 3.05 fim C 2 H 2 /HCN feature. 

Possible explanations of the discrepancy between the model and 
the measurements are discussed in more detail in Sect. 4.3.5. 
This last model is, however, considered as the best dynamic 
model for RScl, among the ones available, and will be used 
in the comparison to the observed spectro-interferometric 
measurements of the star in the following. 



4.3.4. Self-consistent formation of extended molecular 
structures 

As discussed in Sect. 4.2.2, the N-band visibility data display a 
slope for the largest 60 m baseline in the 8-9 fim spectral range 
likely related to the emission of extended molecular structures 
of C2H2 and HCN. Figure 13 shows the comparison of the vis- 
ibility of the best dynamic model for a given cycle with the 
best-fitting COMARCS+DUSTY model superimposed on the 
MIDI visibility data around 60 m for all the available phases: 
0.23/079/0.97. From this figure, we can clearly see that the 
COMARCS+DUSTY model is too optically thin, showing a flat 
visibility profile in the 8-9 //m spectral range (see Sect. 4.2.2). 
On the contrary, the dynamic model is able to reproduce the 
slope of the MIDI visibilities well through this wavelength 
range. This means that the extension of the dynamic model is 
suitable for revealing the presence of the molecular structures 
of C2H2 and HCN, emitting in that spectral range and located 
above the stellar photosphere. However, the dynamic model is 
not optically thick enough to reproduce the measurements in the 
reddest part of the mid-infrared spectral band. 




Wovelength fam) 



Fig. 13. Comparison of the visibility of the best dynamic model (solid 
line) with the best-fitting COMARCS+DUSTY model (dashed line; see 
Sect. 4.2.1) superimposed on the 60 m MIDI visibility data (error bars) 
at phases 0.23/0.79/0.97 (data sets #1, #7, and #21, respectively). The 
yellow zone corresponds to the region dominated by warm molecular 
layers. 

4.3.5. Comparison of the dynamic modeling results to the 
time-dependent spectro-interferometric data 

In the following, we investigate 100 phases of the model 
randomly distributed over 56 cycles corresponding to a time se- 
quence of 60 years. As there is no possibility to define the model 
cycle corresponding to the observations, we decide to follow a 
statistical approach and to represent all the available cycles for 
a phase range close to the observed one (see Sect. 4.3.2). 

Figure 14 shows the time-dependent cycle-to-cycle averaged 
fluxes of the best dynamic model superimposed on the MIDI 
and ISO/SWS spectrometric data in the 8.5, 11.3, and 12.5 fim 
spectral bands. The figure shows that the model agrees relatively 
well with the time-dependent flux data at 8.5 fim, whereas the 
model is fainter than the measurements in the 1 1.3 and 12.5 fim 
spectral bands. 

In compliance with the observations, the comparison of the 
cycle-to-cycle averaged visibilities of the best dynamic model 
from phase-to-phase, considering the same baselines and phases 
as the observations of RScl (see Sect. 3.4), does not show any 
significant variability of the object in the K and N bands. 

Figure 15 shows the cycle-to-cycle averaged visibility of 
the best dynamic model superimposed on the VINCI data at 
phases 0.17 and 0.23. The model shows good agreement with 
the data implying that the extension of the model is suitable in 
the K-band. 

Figure 16 shows the time-dependent cycle-to-cycle averaged 
visibilities of the best dynamic model superimposed on the 
[14-16m], [29-31 m], and [59-62m] baselines MIDI data in 
the 8.5, 11.3, and 12.5 fim spectral bands. The figure only 
presents measurements showing a deviation from the best-fitting 
COMARCS+DUSTY model (see Fig. 10) smaller than the ob- 
ject variability scattering at 15 and 31m baselines (see Fig. 7). 
This allows us to avoid comparing the spherically symmetric 
model with visibility measurements exhibiting a potential 
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Fig. 14. Time-dependent cycle-to-cycle averaged fluxes of the best dynamic model (solid line) superimposed on the MIDI (squares) and ISO/SWS 
(triangles) spectrometric data in the 8.5 (left), 11.3 (middle), and 12.5 yum (right) spectral bands. Dotted error bars give the dispersion created by 
the cycle-to-cycle variation. 
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Fig. 15. Cycle-to-cycle averaged visibility profiles of the best dynamic 
model (green solid line: phase 0.17; blue dashed line: phase 0.23) su- 
perimposed on the VINCI data at phase 0.17 (green squares) and phase 
0.23 (blue triangles). Dotted error bars give the dispersion created by the 
cycle-to-cycle variation. In the inset, a zoom of the relevant frequency 
region is represented. 



departure of the object from sphericity. As we do not have any 
information on the variability effect at the 60 m baseline, all the 
measurements are represented. From this figure, it is obvious 
that the model is not extended enough to reproduce the level 
of the mid-infrared interferometric data. Only the visibility 
measurements at the highest spatial frequencies (B=60m and 
/l=8.0-9.5 /mi), probing regions corresponding to the warm 
extended molecular layers (from 1.5 to 2 R*), are in good 
agreement with the model. 

Besides the possibility of an inaccurate estimation of the pa- 
rameter values of the hydrostatic initial structure, the most prob- 
able reasons for the discrepancy between the dynamic model 
and the spectro-interferometric measurements are the following 
ones. 

(i) The sparse sampling of the transition region from 
windless models to models with considerable outflows in 
certain critical stellar parameters: 

The 100 K higher value of the dynamic model effective tem- 



perature compared to the hydrostatic model could be a reason 
for the discrepancy of the model to the MIDI interferometric 
measurements, given the closeness of the model to the mass-loss 
threshold defined by Mattsson et al. (2010). Another critical 
parameter for dynamic models, in particular for relatively weak 
winds, is the C/O ratio. After testing a model having a C/O 
ratio of 1.47, while all the other parameters are identical to 
the best dynamic model, we found that the model appears, in 
that case, too much extended compared to the 15 m baseline 
VINCI visibility data, whereas the levels of the 60, 30, and 
15 m baselines MIDI visibility data are reproduced relatively 
well. The transition from windless models to models with 
considerable mass-loss rates occurs in a very narrow range of 
stellar parameters especially for the pulsation amplitude, the 
C/O ratio, and the effective temperature, as already pointed out 
by Gautschy-Loidl et al. (2004). In summary, a better sampling 
of parameter space beyond the models currently available 
(Mattsson et al., 2008, 2010), especially in effective temperature 
and C/O, could help reduce the discrepancies between models 
and observations. Such a study is, however, beyond the scope of 
this paper. 

(ii) The approximation of the stellar pulsation with the 
sinusoidal piston: 

As the star is a semi-regular variable showing two main periods 
(374 and 1804 days) with harmonics in the smallest period 
(Whitelock et al., 1997), using the simple sinusoidal piston to 
describe the stellar pulsation could be a crude approximation 
for such a star. The overall wind properties are mostly affected 
by the pulsation amplitude and not so much by the actual form 
of the piston motion (see models based on different types of 
pulsation descriptions presented by Freytag & Hofner 2008). 
However, the particular phase coupling between dynamic and 
radiative properties of the models imposed by the sinusoidal 
piston might cause some of the discrepancies between models 
and observations encountered here. 

(iii) The approximation of a complete momentum and po- 
sition coupling of gas and dust (i.e., absence of drift) and the 
approximation of the small-particle limit (or Rayleigh limit) 
in the determination of the dust grain opacity in the dynamic 
computations: 

Gas and dust decouple when the gas of the wind is diluted as is 
the case for low mass-loss rate stars like R Scl. Sandin & Hofner 
(2003, 2004) find that decoupling the dust and gas phase in- 
creases the dust formation, which could reconcile the dynamic 
model with the mid-infrared spectro-interferometric measure- 
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Fig. 16. Time-dependent cycle-to-cycle averaged visibilities of the best dynamic model (solid line) superimposed on the [14-16 m], [29-31 m], and 
[59-62 m] baselines MIDI data (squares) in the 8.5 (left), 11.3 (middle), and 12.5 yum (right) spectral bands. Dotted error bars give the dispersion 
created by the cycle-to-cycle variation. 



merits. Similar arguments, regarding the difficulty of produc- 
ing low mass-loss rate models, have also been pointed out by 
Mattsson et al. (2010). 

It is found by Mattsson et al. (2010) that dust grains form- 
ing in slow wind atmospheres may grow beyond sizes where 
their opacities can be approximated by the small-particle limit. 
Moreover, it has been demonstrated by Hofner (2008) that ra- 
diation pressure on micron-sized silicate grains may play a key 
role in driving outflows in M-type AGB stars. Similar arguments 
could also apply to low mass-loss rate carbon-rich stars like 
RScl. 



5. Conclusions and perspectives 

In this work, we have presented the first interpretation of 
combined photometric, spectrometric, and interferometric mea- 
surements of a carbon-rich star based on state-of-the-art, self- 
consistent dynamic atmospheric models. 

Interferometric observations do not show any significant vari- 
ability effect at the 16 m baseline between phases 0.17 and 0.23 
in the K band. This is, however, not surprising given the short 
time interval and because the VINCI broadband measurements 
average out the short term effects of a shock front in CO or 
CN lines. No significant variability effect is found for both 15 m 



baseline between phases 0.66 and 0.97, and 31m baselines be- 
tween phases 0.90 and 0.97 in the N band. This means that the 
stellar radiative pressure is not strong enough to reveal a move- 
ment of the warm mass dust shells larger than 3 AU, in good 
agreement with the amplitude of the mass shell of theoretical 
carbon-rich dynamic models. 

The spectro-interferometric predictions of the dynamic model 
atmosphere are in relatively good agreement with the dynamic 
picture that we deduced from observations. We find rather good 
agreement between the dynamic model and the spectrophoto- 
metric data from 0.4 to 25 jim. The model agrees relatively 
well with the time-dependent flux data at 8.5 //m, whereas the 
model is too faint compared to the measurements in the 1 1 .3 and 
12.5 /mi spectral bands. The 15 m baseline VINCI visibilities are 
reproduced well for the two post-maximum brightness phases, 
meaning that the extension of the dynamic model is suitable in 
the K-band. In the mid-infrared, the slope of the 60 m baseline 
MIDI visibilities in the 8-9 //m spectral range is reproduced well 
for all the available phases. This means that the dynamic model 
has the proper extension to reveal the molecular structures of 
C2H2 and HCN located above the stellar photosphere, whereas 
the non consistent COMARCS+DUSTY model fails to do so. 
The discrepancy between the dynamic model and the spectro- 
interferometric data could be related to inaccurate estimation 
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of parameter values of the hydrostatic initial structure, a differ- 
ence between the parameter values deduced from the hydrostatic 
modeling and the ones related to the hydrostatic initial structure 
of the dynamic model, the sinusoidal piston simulating the stel- 
lar pulsation, the inclusion of SiC opacities in the a posteriori 
radiative transfer, or a combination of all these effects. However, 
it seems that, owing to the strong nonequilibrium process of dust 
formation in AGB stars, the transition from models without wind 
to models with considerable mass-loss rates occurs in a very 
narrow parameter range. The most sensitive parameters allow- 
ing this sharp and strong transition between models with and 
without wind are the effective temperature, the amplitude of the 
piston velocity, and the C/O ratio. Therefore, it seems necessary 
to improve the sampling of critical regions in parameter space 
in the grid of hydrodynamic models for further investigations of 
the extended structures of low mass-loss carbon stars like R Scl. 
Finally, the complete dynamic coupling of gas and dust, and 
the approximation of grain opacities with the small-particle 
limit in the dynamic calculation, both probably unsuitable for 
low mass-loss rate object like RScl, could also contribute to 
the disagreement between the dynamic model and the spectro- 
interferometric data. First tests with C-star models based on 
opacities that take grain size into consideration (Mattsson & 
Hofner, in prep.) show that wind characteristics may be affected 
considerably in models close to a mass-loss threshold (Mattsson 
et al., 2010), such as the one used here to represent the low 
mass-loss C-star R Scl, resembling recent wind models for M- 
type AGB stars (Hofner, 2008). 
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